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ABSTRACT 

Context. Unbiased molecular line surveys are a powerful tool for analyzing the physical and chemical parameters of astronomical objects and 
are the only means for obtaining a complete view of the molecular inventory for a given source. The present work stands for the first such 
investigation of a photon-dominated region. 

Aims. The first results of an ongoing millimeter-wave survey obtained towards the Orion Bar are reported. 

Methods. The APEX telescope in combination with the APEX-2A facility receiver was employed in this investigation. 

Results. We derived the physical parameters of the gas through LVG analyses of the methanol and formaldehyde data. Information on the 
sulfur and deuterium chemistry of photon-dominated regions is obtained from detections of several sulfur-bearing molecules and DCN. 
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1. Introduction 

Photon-dominated regions (PDRs) are commonly found inter- 
stellar environments whose global properties are determined by 
intense far-ultraviol et (FUV) radiation em erging from nearby 
young OB stars (e.g. Holle nbach & Tielensll997l) . Owing to its 
proximity and nearly edge-on orientation, the Orion Bar PDR 
has received particular attention as a template for studies of the 
spatial stratification of various atomic and molecular species 
from the highly penetrated surface layers deep into the parental 
molecular cloud. Targeted (sub)millimeter-wave investigations 
indicate that, besides a considerable number of ubiquitous as- 
tronomical molecules such as CO, CS, HCN, CH3OH, and 
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pears to be CF + , which so far has only been detected toward 
the Orion Bar (Neufeld et al 2006, this volume). PDRs are 
also thought to show enhanced abundances of molecules car- 
rying refractory ele ments due to grain breakup caused by the 
FUV radiation field dSchilke et allkOOll) . Additionally, obser- 
vations of the Horsehead nebula suggest that at least part of the 
molecular carbon chain budget in PDRs may be prod uced from 
UV-destruction of PAHs and carbonaceous grains dPetv et alJ 
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In view of this, PDRs are attractive targets for unbiased 
molecular line studies, which also help to derive the global pic- 
ture of the physics and chemistry associated with them. Here, 
we present the initial results from an unbiased molecular line 
survey obtained between 279 and 308 GHz with the APEX tele- 
scope (Giisten et al. 2006, this volume). 

2. Observation 

Observations were conducted with the APEX 12-m telescope 1 
towards the Orion Bar at the position a2ooo = 5 /! 35 m 25.3 s , 
£2000 = -5 °24'34.0", correspon ding to the "Orion Bar (HCN)" 
position of Schilke'et alJ (|200ll). the most massive clump seen 
in H 13 CN JLis & Schilkell2003l hereafter LS03). Data were 
taken in 2005, between October and December in the position- 
switching mode, with the reference position at (600", 0"). The 
APEX-2A receiver (Risacher et al. 2006, this volume) was used 
in combination with the Fast Fourier Transform Spectrometer 
(Klein et al. 2006, this volume) providing a bandwidth of 
1 GHz and a resolution of 0. 12 km s _1 . The frequency range be- 
tween 279 and 307.7 GHz was completely covered; additional 
data were taken in selected frequency ranges from 318.5 GHz 
to 361.5 GHz. Step- widths of approximately 500 MHz were 
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1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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used to cover each frequency setup twice to facilitate the side- 
band assignment. The antenna temperature was converted to 
the main-beam brightness temperature by using forward and 
beam efficiencies of 0.97 and 0.74, respectively. The noise 
level is of the order of 0.06 K-0. 1 K, and the system temper- 
ature around 160 K. The calibration was performed by using 
the APECS software (Muders et al. 2006, this volume). The 
pointing was checked on Mars and was found to be accurate 
within a few arc-seconds. We estimate a calibration uncertainty 
of ~ 30%. 

3. Analysis of the data 

The line identification was b ased on the Colo gne Database 
for M olecular Spectrosco py 2 (LVIuller et alJ 120011) and on the 
JPL 3 JPickett et alJ Il998l) line catalog and performed with 
the XCLASS software. We identified 16 different species and 
a number of their isotopologs; no strong unidentified lines 
were found in the range surveyed. In the following para- 
graphs, we focus on the excitation conditions of methanol 
(CH3OH) and formaldehyde (H2CO), on the one hand, and 
on the analysis of S-bearing molecules, on the other. Methanol 
and formaldehyde are present in our dataset with a large num- 
ber of tr ansition s (Table I D, and their c ollisional rates are 
available (Gre enlll99lUPottage et alJ l2002). Thus, a more rig- 
orous multi-line analysis is possible for these two species, 
which delivers temperature and density estimates. As for the 
S-bearing molecules in our dataset, some of them, namely 
CS, H 2 S, SO, SO +, and SO?, have already been studied to- 
wards this source dHoeerheiide et al1ll995t Ijansen et al.ll 19951: 
iFuente et al.ll2003h . whereas HCS + , H 2 CS and NS are being 
analyzed here for the first time in this environment. A compre- 
hensive analysis will be presented once the line survey has been 
completed. 

3.1. CH 3 OH and H 2 CO 

The excitation analysis of CH3 OH and H?CO was p erformed 
by using the method described bv lLeurini et alJj2004l) . recently 
modified to treat H2CO as well. The procedure is based on 
the simultaneous fit of all the lines with a synthetic spectrum 
computed in the LVG approximation with the cosmic back- 
ground as the only external radiation field. If lines are optically 
thick, the model fits source size, temperature, H2 density, and 
CH3OH/H2CO column density; for optically thin lines, source 
size and column density cannot be determined independently, 
and beam-averaged column densities are provided. 

To improve the signal-to-noise ratio in the weaker transi- 
tions, we smoothed the data to a resolution of 0.5 km s _1 . At 
this resolution, each line has a single-peak Gaussian profile 
(v lsr = 10.0 ± 0.2 km s -1 , Av = 1.7 ± 0.3 km s" 1 ), while, at 
the original resolution of 0.12 km s _1 , two peaks were detected 
in the 5i, 5 -> 4 lA H 2 CO and in the 6-1 -» 5_i CH 3 OH-£ 
transitions. Non-Gaussian profiles were found in other transi- 
tions as well. By fitting the 6_i — > 5_i CH30H-£ line with 

2 http://www.cdms.de 

3 http://spec.jpl.nasa.gov 



Table 1. Line parameters of the observed transitions 



transition 


frequency 




J T MB dv 




[MHz] 


[K] 


[Kkms- 1 ] 


H 2 C0 4 li4 -> 3 U 


281526.929 


32.06 


10.7 


±0.5 


H 2 CO 4 , 4 -» 3 ,3 


290623.405 


20.96 


5.9 


±0.4 


H 2 CO 4 2 3 -> 3 2 ,2 


291237.767 


68.09 


2.4 


±0.3 


H 2 CO 4 3>2 -» 3 3 ,i 


291380.488 


126.95 


3.4 


±0.3 


H 2 CO 4 3>1 -> 3 3 , 


291384.264 


126.95 


3.7 


±0.3 


H 2 C04 U -> 3 1-2 


300836.635 


33.45 


9.8 


±0.5 


H 2 C0 5^ 5 -» 4i, 4 


351768.645 


45.57 


10.1 


±0.9 


CH 3 OH-£ 6 -> 5 


289939.477 


47.93 


0.5 


±0.1 


CH 3 OH-£6_! -> 5_j 


290069.824 


40.39 


0.9 


±0.2 


CH 3 OH-A 6 -> 5 


290110.666 


34.82 


1.2 


±0.2 


CH 3 OH-£6! -> 5i 


290248.762 


55.87 


0.3 


±0.2 


CH 3 OH-£ 6 2 -» 5 2 


290307.376 


60.72 


0.3" 


±0.2 


CH 3 OH-£ 6-2 -> 5-2 


290307.643 


57.07 






CH 3 OH-A U -> h 


303366.890 


2.32 


0.6 


±0.2 


CH 3 OH-A 2 1 -> 2 


304208.350 


6.96 


0.8 


±0.2 


CH 3 OH-A3i -> 3 


305473.520 


13.93 


0.9 


±0.2 


CH 3 OH-A 4 1 -» 4 


307165.940 


23.21 


1.4 


±0.3 


C 34 S 6-^5 


289209.068 


34.70 


4.6 


±0.7 


C 33 S 6-^5 


291485.935 


34.98 


1.6 


±0.3 


CS6^ 5 


293912.087 


35.27 


40.7 


±2.2 


CS7^6 


342882.850 


49.37 


33.1 


± 1.9 


HCS + 7^6 


298690.453 


43.01 


0.3 


±0.1 


HCS + 8^7 


341350.229 


57.34 


0.6 


±0.1 


SO 7 6 -> 65 


296550.064 


50.66 


3.7 


±0.5 


SO 7 7 -> 6 6 


301286.124 


56.50 


4.0 


±0.6 


SO 7 8 -> 6 7 


304077.844 


47.55 


7.2 


± 1.3 


SO 87 -> 7 6 


340714.155 


64.89 


3.6 


±0.7 


SO 8 8 -> 7 7 


344310.612 


70.96 


3.0 


± 1.1 


SO 89 -> 7 8 


346528.481 


62.14 


7.2 


± 1.2 


5fi+ 2tt 1 n /<-) . 1 1 try 

oU Hi/2 1-3/.Z — * 1 1 / ze 


Qni QAi 
JUL 30L.JUL 


Jo.yZ 


0.6 


±0.1 


so + 2 n 1/2 13/2 -> 11/2/ 


301736.791 


39.03 


0.4 


±0.1 


H2CS 9i 9 — * 8i s 


304307.645 


71.60 


0.6 


±0.1 


H2CS 10o.l,o — > %,9 


342946.335 


74.13 


0.4 


±0.1 


H 2 CS 102,9 -> 9 2 8 


343322.111 


126.83 


0.4 


±0.1 


H 2 CS IO1.9 -> 9 18 


348534.225 


88.47 


0.8 


±0.2 


NS 2 n 1/2 13/2 -» ll/2e 


299700.097 


39.82 


0.4 


±0.1 


NS 2 n 1/2 13/2 11/2/ 


300098.057 


39.93 


0.4± 0.1 


NS 2 n 1/2 15/2 -> 13/2e 


345823.532 


54.20 


0.7 


±0.2 


NS 2 n 1/2 15/2 -> 13/2/ 


346220.774 


54.33 


0.8 


±0.2 


H2S 3 3 ,o — ^ 32,i 


300505.560 


154.48 


0.5 


±0.1 



" blended with the 6_ 2 -> 5_ 2 CH3OH--E line. 



the GAUSS method of CLASS, we found two velocity com- 
ponents at visr = 10.89 ± 0.06 and v Lsr = 9.82 + 0.05 km s~', 
with Av = 0.62 + 0.13 and Av = 0.63 + 0.06 km s _1 . These 
values resemble what LS03 found for the H 13 CN clumps 1, 5 
and 8, which are all within our beam. For the sake of simplicity, 
one velocity component was used both for the CH3OH and the 
H2CO analysis, as derived from the smoothed CH3OH data. 

High density (« ~ 7 x 10 6 cm 4 ) and a moderately high 
temperature (T ~ 70 K) are needed to reproduce the CH3OH 
spectrum. The temperature is constrained to this value by the 
non-detection of high-excitation (E/ ou . > 80 K) lines (Fig.0. 
The high density value we derived suggests that CH3OH comes 
from small clumps. Therefore, guided by LS03, we used a 
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Fig. 1. Spectrum of the 6k — > 5,t CH3OH band. Overlaid on the 
data, in blue, is the synthetic spectrum corresponding to our 
best fit. 

source size of 5" and derived N(CH 3 OH-£) = 7 x 10 14 ctrT 2 
and N(CH 3 OH-A) = 8 x 10 14 cm" 2 . Our estimate of n in 
the c lumps agrees with previous studies (e.g. lYoung Owl et alJ 
2000; LS03); Tt,-„ is also similar to previous results for the 
clu mp gas (see discussion in LS03) and to theoretical mod- 
els jGorti & HollenbachlEooi) . The H 2 column density (9 x 
10 22 cm 2 ) was derived from our observations of the C 17 Q(3- 
2) line, assuming 17 0/ 16 0~ 1790 JWilson & Roodll994 and 
T ro t=7Q K. The methanol column density then corresponds to 
a fractional abundance averaged over the beam of 10~ 9 , typical 
of dark clouds. The 3<x confidence level for Tki„ ranges between 
50 and 75 K, and (5-20)x 10 6 cirT 3 for «h 2 . 

The parameters derived from H2CO differ significantly 
from the ones of CH3OH, as high temperature (T ~ 150 K) 
and moderate density (« ~ 5 x 10 5 cm 4 ) are needed to re- 
produce the data. Since all H2CO lines are optically thin, we 
derived a beam-averaged column density of 9 x 10 13 cirT 2 for 
ortho-H2CO and 3 x 10 13 cirT 2 for para-H?C O, correspond- 
ing to a fractional abundance of 10~ 9 . iFuente eta D J! 994 
found the same value, although at another position in the 
Bar. Our analysis infers a Tun in agreement with the re- 
sults of lBatrla & WilsorJ J2003I) from NH^ b ut higher than 
the one derived bv iHogerheiide et alJ dl995h . 85 + 30 K, 
also with formaldehyde, but toward a different position. The 
density derived is substantially higher than that proposed by 
IHogerheiide et alJ i 19951) for the interclump m edium (3 x 
10 4 cm " 3 ), but it agrees with the value inferred bv lSimon et alJ 
ll 1997b based on CN. While the «h 2 is well-constrained (1— 
5xl0 5 cirT 3 ), the temperature we derived is a lower limit. 

The infrared pumping through an internal radiation field 
would probably affect our results for both molecular species. 
However, since no model of the IR field is available for the 
Bar, no internal radiation field was taken into account. 

Our data suggest that CH3OH is found mainly in the 
clumpy medium, while H 2 CO traces the interclump mate- 
rial. Both methanol and formaldehyde can form on grain sur- 
faces (Watanabe & Kouchil 12002). It therefore seems remark- 
able that they trace different environments. We speculate that 
in the interclump gas, methanol, once released from grain 
surfaces, is photodissociated to form formaldehyd e (one o f 
the possible photodissociation products, ILe Teuff et al. 2000). 



Alternatively, our observations may indicate a different for- 
mation mechanism for the two species, surface ch emistry for 
CH3OH (as no gas-reaction can efficiently form it, iLuca et all 
2002), and gas phase chemistry for H2CO (via th e neutral- 
neutral reaction CH 3 + O, e.g. lLe Teuff et aljEo OO). More de- 
tailed modeling, including release of molecules from grain sur- 
faces in the PDR, is clearly needed. 

3.2. Sulfur-bearing species 

The excitation analysis of the S-bearing molecules listed in 
Table \l\ was performed by fitting all the observed transitions 
with a synthetic spectrum computed unde r the LTE a ssump- 
tion, in the manner described bv IComito e t al. (2005). Since 
the kinetic temperature, the source size, and the column den- 
sity are degenerate parameters in the optically thin limit, we 
fixed Tkin to the value derived from H2CO and fit the beam- 
averaged column density of each species. For SO2, OCS, and 
HCS, upper limits to the column densities were derived, based 
on the non-detections of lines in our survey range. Results are 
given in Tabled The uncertainties on the column densities are 
derived with a^- 2 analysis and correspond to the 3cr confidence 
level. CS and its isotopologs are not included in our analysis, 
because of the an omalous 33 S/ 34 S ratio of ~ 3 we observed (for 
comparison, see IChin et alJl 996). Observations of other rota- 
tional transitions are needed to constrain this ratio reliably. 

ISternberg & Dalgarnol ll 19951 hereafter SD95) studied the 
chemistry of PDRs produced in a dense molecular cloud ex- 
posed to intense far-UV radiation fields. The density they em- 
ploy (tih = 10 6 cirT 3 ) is higher than the one relevant here for 
the interclump gas; discrepancies in the results may be due to 
this. Their chemical network does not include NS and H 2 CS, 
but all the other S-bearing molecules observed here are found 
there. 

According to SD95, the SO + /SO ratio (observed 0.2 from 
Table |3 can be used as diagnostic of the different chemical 
zones in the cloud and changes between 0.23 at A v = 1.5 and 
0.0028 at A v = 3. Assuming that the Bar is edge-on, based 
on the distance between the ionization front and our position 
(~ 30"), and assuming an H2 density of 10 5 cm -3 , we de - 
rived a visual extinction of > 10 mag llFrerkingetalJll982l) . 
Given the complex geometry of the Orion Bar, which has edge- 
on (but oblique) and face-on parts (Hogerheiide et alJl 19951) . it 
is likely that together with the molecular material at this ex- 
tinction, a hotter, more diffuse layer of gas, corresponding to a 
face-on surface layer, contributes to the emission we observe. 
SD95 predict S0 2 to be more abundant than SO and SO + at 
high A v , while our upper limit indicates that SO2 is less abun- 
dant. Similar to the result found for SO + /SO, our estimates of 
SO2/SO and S02/SO + suggest that these species have a signif- 
icant contribution from a hotter layer of gas with lower extinc- 
tion than the molecular layer. 

IFuente et alJ J2003I) report different column densities for 
SO + and SO2; however, their SO + /S02 ratio ranges between 
1.2 at the IR front and 0. 4 at (20",-20 " ) from it, which is sim- 
ilar to our value of 0.3. IJansen et alJ i 19951) found values of 
N(H 2 S) between 8.9 x 10 12 cm™ 2 and 3.9 X 10 14 cm™ 2 across 
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Fig. 2. H 13 CN(4-3) and DCN(4-3) lines. For clarity, the H 13 CN 
line has been translated by 0.2 along the y-axis. 

the Bar, with our value of 1.3 x 10 13 ctrT 2 lying in between. 
Discrepancies in the absolute values of the column densities 
may depend on the source sizes, on temperatures, and on the 
total H2 column density at the observed position. 



Table 2. Column densities of S -bearing molecules 



Species 


N [cnT 2 ] 


Species 


N [cm" 2 ] 


NS 


(1.3 ±0.3)10" 


so 2 


< 3 10 13 


SO 


(1.0±0.2)10 14 


HCS 


<4 10 14 


so + 


(2.0±0.5)10 13 


OCS 


< 4 10 13 


H 2 S 


(2.5 ± 1.0)10 13 






H 2 CS 


(2.5 ±0.4) 10 13 






HCS + 


(2.5 ±0.7)10 12 







4. Detection of a deuterated molecule : DCN 

Molecular deuteration studies have been extremely powerful 
in probing the physical history of sources. The deuterium 
isotopic ratio in molecules can indeed be enhanced in low- 
temperat ure environme nts compared to the cosmic D/H ratio 
(~1Q- 5 . lrlnskvl 1200^1 . This enhancement proceeds initially 
from the transfer of deuterium from the main reservoir HD to 
two reactive ions (H2D + a nd CH?D + ) through exo thermic re- 
actions with H+ and CH 3 + (Roberts & Mill arl20od and refer- 
ences therein). The endothermicity of the backward reactions 
(respectively by 220 K and 370 K) ensures that the H 2 D + /H+ 
and CH2D + /CH^ ratios are significantly enhanced at low tem- 
peratures. Deuterium is then channeled to other molecules by 
gas-phase reactions with these two ions and by grain surface re- 
actions. Highly-deuterated molecules can also be found in hot 
gas, where they are then out-of-equilibrium fossils of an ear- 
lier cold phase (e.g. in hot corinos around low-mass protostars, 
where they are believe d to be remnants of the p restellar cold 
and depleted chemistrv: IParise et al 120041 12OO6). 

The DCN(4-3) transition at 289.645 GHz (E up = 34.56 K) 
was detected in our survey, and found to be very bright (Fig.|4j. 
This detection, together with the detection of DCO + in the 
Horsehead nebula (Pety et al. 2006, in preparation), is opening 
up the study of fractionation processes in PDRs. The H 13 CN(4- 
3) transition (E up = 4 1 .76 K) was also detected at 345 . 340 GHz. 
The H 13 CN emission was shown from interferometric observa- 



tions (LS03) to be mostly associated with clumps in the PDR. 
The DCN emission is thus very likely to also be associated with 
clumps, as DCN would not survive the high temperatures of the 
interclump gas. The linewidths of DCN(4-3) and H 13 CN(4-3) 
are comparable and non-Gaussian profiles are detected in both 
cases. We thus assume that H 13 CN and DCN trace the same 
region. As the energies of the two lines are very similar, the de- 
rived abundance ratio is nearly independent of the temperature. 
Assuming LTE and a 12 C/ 13 C ratio of 70, we derived beam- 
averaged DCN and H 13 CN column densities between 8.4x 10 11 
and 1 .7x 10 12 cm" 2 , and between 1 .7x 10 12 and 2.8x 10 12 cm" 2 , 
respectively, in the range of temperatures 30-150 K. This cor- 
responds to a DCN/HCN ratio of 0.7-0.9%, which is an inter- 
mediate val ue between that obs erved in warm gas (Orion Hot 
Core, 0.1%. ISchilkeet alJI 992h on the one hand, and in dark 
clouds (LI 34. 5%. lTurnerl200lh or cold (30- 50K) gas of the 
OMC 1 ridge region (1-6% ISchilke etal.l l 992). on the other. 

In hot cores/corinos, where the gas has not yet had time to 
return to equilibrium (which would require timescales of ~ 10 4 
yrs), abundant deuterated molecules can be considered as rem- 
nants of earlier cold prestellar chemistry, most probably stored 
in the ices on grain surfaces and then released when the proto- 
star heats its surroundings. In contrast, the gas in the shielded 
clumps of the bar might be in steady state. In this case, the 
observed deuterated molecules have to be present-day product 
molecules. The detection of DCN in this relatively warm envi- 
ronment, where H 2 D + would not survive, is in this sense quite 
striking and is a strong hint of the efficiency of the fractionation 
processes occurring via the CH 2 D + channel, which survives 
higher temperatures than H2D 4 ", due to the h igher exothermic- 
ity of the CH^ reaction with HD (Turner 2001). The Orion Bar 
thus appears a unique reference to test the fractionation reac- 
tions involving CH2D + . Knowing the spatial distribution of the 
DCN/HCN isotopic ratio would then be an interesting tool for 
understanding the effect of extinction on the fractionation pro- 
cesses. 

Acknowledgements. BP is grateful to the Alexander von Humboldt 
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